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ABSTRACT 



A century of observations of u Centauri allows us to accurately follow the period 
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changes of a large sample of variable stars. Although the period changes are composed 
of monotone and irregular changes, period increases dominate in most of the RRab 
stars. The typical period increase rates are of the same order of magnitude as model 
calculations predict for the evolved, redward phase of horizontal branch evolution. For 
• the 44 well observed RRab stars which show monotonic period changes, a mean rate of 

0.15 d/Myr has been found. The period changes of the first overtone RRc type stars 
show a much more complex, irregular behavior. 
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1. Introduction 
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One of the principal goals of globular cluster photometric studies carried out during the past 
century was to determine period changes of RR Lyrae stars. The hope to reveal evolutionary 
changes motivated many extensive works. (For details see e.g. Smith (1995) and Rathbun and 
Smith (1997).) However, attempts to connect the observed period changes directly with the different 
phases of horizontal branch (HB) evolution have not led to indisputable conclusions. Neither the 
large period decreases nor the random and/or abrupt period changes that have been detected in 
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many clusters can be explained on the basis of steady evolutionary effects. 

The basic tool in the determination of period changes of variable stars is the O — C diagram. 
This plots the difference between the observed times (O) of a particular phase (e.g. maximum) of 
the light-curve and the predicted times of the same phase (C) calculated according to an accepted 
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ephemeris. A nonlinear O — C indicates period changes while a linear O — C with a non-zero slope 
points to inaccuracy of the adopted period. 

Sterne (1934) showed that period-noise can have a cumulative effect in the O — C of variable 
stars, thus mimicking true period changes. However, in case of RR Lyrae stars Balazs-Detre and 
Detre (1965) concluded that their O — C diagrams cannot be wholly explained as the result of 
period-noise. Nevertheless, Lee (1991) showed that the mean period change rates observed in 
globular clusters were consistent with synthetic HB models if the HB types of the clusters were 
also considered. However, accurate determinations of period changes are available for only a small 
sample of clusters and for many of these, the derived period change rates vary from one study to 
another (see e. g. the 0.07 d/Myr difference in the mean period change rate for NGC 7006 given by 
Wehlau et al. (1992) and Wehlau et al. (1999)). Thus it is still an important task to define accurate 
period change rates in globular clusters by using not only the most recent data, but also any 
earlier observations, even if they can only be reconstructed from published light curves or epochs 
of maximum light. Searching for unpublished cluster observations and making them accessible for 
future studies is very worthwhile. 

The preference for period increases in u> Centauri was first noticed by Martin (1938). Belserene 
(1964) strengthened this result and raised the possibility of an evolutionary explanation. The recent 
CCD observations of the variables (Kaluzny et al. 1997) within the framework of the OGLE (Optical 
Gravitational Lensing Experiment) project have lengthened the time-base of the observations to 
more than one hundred years, thus making it possible to extend the earlier period change studies. 

2. Data 

We used all known and available photometric information on u Centauri variables including the 
original photographic data of Bailey (1902); Martin (1938); Dickens and Saunders (1965); Dickens 
and Carey (1967) and the recent CCD V observations of Kaluzny et al. (1997). In addition, the 
scarce photoelectric V observations made by Sturch (1976, 1978) and Eggen (1961) of V92 proved 
to be useful in some cases for defining anO-C value at an epoch when no other observations were 
obtained. 

Belserene (1964) published O — C values for 47 RRab stars based partly on observations 
obtained with the Yale-Columbia 26-inch refractor in 1947, 1948, 1956, 1958, 1961, 1962 and 1963. 
The individual magnitudes were not included in her paper. However, we have obtained most of 
the unpublished data for these 47 RRab stars and for a further 41 variables, mainly RRc stars. In 
addition, magnitudes derived from observations made with the same telescope in 1965 and 1967 
have been obtained. In Table 1, a sample of the magnitudes for the first 11 stars are listed; the rest 
are available in the electronic edition. Observations of some variables not included in the present 
study are also given in Table 1 (SX Phe: V65, red variable: V164, and variables with not enough 
data to investigate period changes: V167 and V169). We used O — C values that Belserene (1964) 
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derived from the 1922 — 6 observations (Wright 1940) and for the Yale-Columbia 1961 and 1962 
observing seasons because these original data could not be located. 

Wehlau and Bohlender (1982) examined the period changes of BL Her type stars in globular 
clusters. To give the most complete sample of the period changes, we have also included the 
u Centauri stars they studied in our present work, even though new data have not been obtained 
for all of the variables. 

About 400 photographic B plates were taken by C. Clement with the University of Toronto's 61 
cm (f/ 15) Helen Sawyer Hogg telescope at the Las Campanas Observatory of the Carnegie Institu- 
tion of Washington from 1972 to 1986. The plates were measured on a Cuffey iris astrophotometer. 
As these observations span more than a decade, they have crucial importance for following the 
period changes in the second half of the century. In Table 2, part of these data for the first 14 
variables are given. The entire data set is presented in the electronic edition. 

Photographic B and V observations obtained by E. H. Geyer with the 80/90/300cm ADH 
Baker Schmidt camera of the Boyden Observatory in the 1962 observational season (Geyer and 
Szeidl 1970) were also used. The photometric data of these observations have not been published, 
but recovery of the original measurements is possible for 18 variables out of which 16 are used in 
this work. The other two stars are V65 (SX Phe), and V167. Table 3 lists these photographic B 
data. For most of the observed stars only folded light-curves are available. These folded curves 
could be used to determine phase relation only if the assumed period and initial epoch were known 
without doubt. 

In Table 4 the different observations of the u Centauri variables are summarized. The new 
variables discovered by Kaluzny et al. (1997) are not included in Table 4 because there is no previous 
record for determining the period variation. RR Lyraes with not enough data to construct any reli- 
able segment of their O-C are also omitted (V80, 110, 114,135, 136, 137,141, 143, 145, 146, 147,153,154, 
156,157,158,158,169 and all the variables after V170). 

3. Phase shift diagrams, period changes and folded light-curves 

We used the original photometric data to determine phase shifts whenever they were available. 
(These sources are denoted by X in Table 4.) The CCD V (Kaluzny et al. 1997) observations were 
used to define the normal curves and the phase shifts for the other epochs were measured relative 
to these curves. Both vertical and horizontal shifts were derived for all the other photometric 
data in order to match the normal curve best. The most deviant magnitudes were omitted from 
each dataset. We applied a single vertical (zero-point) shift to the different observations, with the 
exception of the OGLE measurements. When the same star appeared in different OGLE fields, 
these data were treated separately and differences in their zero points were also determined. 

For determining O — C values (horizontal shifts) the datasets were divided into shorter segments 
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(1 — 3 years long). (For simplicity, we denote O — C as, the phase shift between the light-curves as 
a whole instead of the traditional definition of reading phase shifts only between special parts of 
the light-curves.) If there were no OGLE observations of the star, the phase shifts were measured 
relative to Martin's (1938) photographic m pg curves. 

As phase shifts between photographic B,V, photoelectric and CCD V observations were mea- 
sured, it was also important to check the effect of data inhomogenity on the results. By determining 
the optimal horizontal and vertical shifts between the B and V light-curves of the recent CCD B 
and V observations of the M2 variables (Lee and Carney 1999), a mean phase shift of 0.0021 d 
(with 0.0039 rms) between the V and B light-curves of the 19 RRab stars was found. The shift for 
the 12 RRc stars is 0.0001 d (0.0021 rms). Consequently, although the times of maxima system- 
atically deviate in B and V, when the phase shifts of the whole light-curves are considered then 
no significant difference according to the colors occur. The O — C ambiguity arising from data 
inhomogeneity is therefore less than 1 — 2% of the period. 

O — C values were derived from the data compiled in Table 4 at all of the epochs for which 
it was possible to obtain a reliable value. Observations of Bailey (1902) and Martin (1938) were 
divided into two parts, while the Las Campanas observations allowed for the determination of six 
independent O — C points in most cases. Observations from the different OGLE fields were used 
independently. If the observations of one of the OGLE fields were extensive enough, O — C -s were 
determined for two of its segments. The O — C values given by Belserene (1964) for 1922 — 6, 1961 
and 1962 (JD 2423530, 2437454 and 2437727, respectively) were appropriately transformed if the 
O — C curve was constructed with a period different from the one she used. Similar transformations 
were applied for the ADH Schmidt data in cases where only the folded light-curves were available. 

We constructed the O — C curve by using an estimated value of the average period over the 
century (P a ) in order to obtain a symmetrical form and fitted the O — C points by different order 
polynomials. The O — C can be expressed as 

0-C = j P{E)dE - EP a = -^j P ^ dt ~ f W 
where E is the epoch number, and P a is the period used to determine O — C . 

If we approximate the O — C with a polynomial in the form of 

k 

0-C = J>f-\ (2) 

i=i 

then the temporal period can be determined using the q coefficients: 

k 

P(t) = P a ~ l)cit l ~ 2 + Pa- (3) 

If it is possible to determine the period with sufficient accuracy at the date of the observations, we 
also have a direct measure of the period changes for comparing with the above fit. 
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Furthermore, if we transform the time according to 

k 

t' = t-Y,Cit-\ (4) 

i=i 

then all the data can be folded with a single period because, from Eq. 2, it follows that the O — C 
term vanishes. 

Thus we have two methods for checking the validity of both the O — C values and their trial 
polynomial fits, namely (a) a comparison of the direct period determinations with the temporal 
periods predicted according to the derivative of the O — C and (b) the scatter of the folded curve of 
the time transformed data. This is essential when an ambiguity in the cycle number counts occurs. 

Figs, la-1, 2a-h and 3 show the results for RRab, RRc and RR Lyrae-like stars (which were 
formerly called BL Her type variables), respectively. In the left panels which show the O — C data 
and their polynomial fits, the filled circles denote O — C points derived from the phase shifts of 
the original data while open circles are for Belserenes's (1964) O — C values for 1922 — 6, 1961 and 
1962 and for phase shifts which were read from Geyer's folded light-curves. If the above checks 
could not help in deciding cycle numbers, and/or there remained significant uncertainties in the 
O — C curve or in its polynomial fit, the fit is shown by a dashed line. The O — C were divided into 
two overlapping parts which were fitted separately with different order polynomials in the cases of 
some RRc stars showing very complex O — C (V10, 64, 76 ,95, 126). This procedure avoids strong 
discontinuities in the O — C and period fits. However, the possibility that the period changes are 
in reality abrupt cannot be excluded. 

In the middle panels, period changes calculated from the O — C fit according to Eq. 3 (dashed 
lines) are shown and compared with the results of direct period determinations. Whenever the 
observations span over two or more consecutive seasons, then direct period determinations with 
sufficient accuracy may also become possible. In some cases periods were determined from single 
observing seasons, but usually the accuracy of these period determinations was not good enough 
to make any use of them. When no ambiguity in the O — C values has arisen, i.e. P a seemed to be 
correct, this period was taken as an initial trial period for each data set. Periods in the vicinity of P a 
were checked and the one that gave the light-curve with the smallest rms scatter in a least squares 
solution was adopted. In those cases when P a might be questioned, i.e. cycle number ambiguity 
has arisen indicating that aliasing might seriously affect the period determinations, a period search 
using the MUFRAN package (Kollath 1990) was performed. If there was more than one possible 
peak in the amplitude spectra of the Discrete Fourier Transform (DFT), the peak which was closest 
to the period determined from the other observations was chosen if it was consistent with the O — C 
fit as well. 

The error range of the periods was determined by checking the residuals (rms) of the folded 
curves using periods close to the one accepted according to the procedure described above. The 
upper and lower limits were set to the periods for which the scatter of the folded curves increased 
by 10%. For the best light-curves (some of Martin's and Kaluzny's observations if more than 200 
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data-points were available, and the quality of the light-curve was high enough to fit a higher order 
(> 7) Fourier sum), only 5% was allowed. These percentages can be justified as the 90% confidence 
region of the period assuming normal errors, correspond e.g. to 5%, 8% and 10% increases in the 
rms scatter of the fit in cases of N=250, M=15; N=120, M=10, and N=60, M=5 (N:number of 
data, M: order of the Fourier fit), respectively (see details in Press et al. (1992)). 

Checking whether the periods determined from Eq. 3 matched the actually observed periods 
within their error ranges could help in many cases to decide the proper cycle numbers. 

The folded curves of all the photometric data derived from the time-transformation of Eq. 4 
are shown in the right panels of Figs. 1-3. As a result of this transformation all the data could 
be brought into phase with the normal curve. In those cases where the cycle number or the shape 
of the O — C fit was dubious, the solution that was accepted was the one that resulted in smaller 
scatter on the folded curve of the time-transformed data. We recall that these folded curves were 
matched to CCD V or to photographic m pg normal curves depending on whether the star was 
measured by the OGLE project or not. 

The scatter of the time transformed data and the agreement between the observed periods and 
the derivative of the O — C fit served as checks on the correctness of the O — C and its polynomial 
fit. Different order polynomials with different weightings, interpolations and extrapolations were 
tested in order to get the most reliable solution. The aim was to determine a more precise O — C 
curve which also enabled us to estimate the contribution of random changes to the global period 
changing term. 

In Table 5 the O — C values, the order of the O — C fit and the actual periods measured and 
their possible ranges at different epochs are given for V3. The results for all of the other variables 
studied are available in the electronic edition which presents Table 5 in its complete form. 

The period change rates defined as f3 =< > are determined by using the polynomial fits 
to the O — C curves. Temporal periods at At = 1000 d intervals are calculated according to Eq. 3 
and the mean of the resultant values is taken as 0. 

The period change rates defined as (3 =< > are calculated by using the polynomial fits 
to the O — C curves. From the temporal periods obtained according to Eq. 3, (At = 1000 d) 
values are determined and their mean is taken as j3. Measuring (5 in this way gives practically the 
same result as calculating 2c^P a from a quadratic O — C fit where C3 is the coefficient of the second 
order term or from measuring the period difference between the first and last observed epochs. 

It is difficult to determine whether a period change is actually random or monotone. However, 
the rms scatter of the calculated values during the century ((tap) may give some indication. 
If the scatter is large, the changes are probably random. In the discussion that follows, we assume 
that no significantly monotonic period change occurs if 

AP 

I < -r— > I < 2a ap i.e. (tap/\/3\ > 0.5. (5) 

At At At 
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P values (in 10 10 days/day and in days/Myr units), the normalized period variations a = j3/ P a 
in 10 -10 /day, (tap, and the <7ap /\(3\ ratio for all the variables studied are given in Table 6. A 

At At 

comparison of these values with Figs. 1-3 demonstrates that our choice of selecting monotone 
period changes according to the Eq. 5 criterion gives reasonable results. 

An examination of the light-curves of all the variables in the course of the above procedure 
revealed that Blazhko behavior (amplitude and/or phase modulation) is manifested in many cases. 
Stars showing definite modulations are marked in Figs. 1-3 and are encoded in column 3 of Table 6. 
Column 4 of Table 6 indicates ambiguous O — C fits, while column 5 refers to the metallicity 
according to Jurcsik (1998b, see also Table 7). 

3.1. RRab stars 

The present investigation confirms Belserene's (1964) result that most of the u> Centauri RRab 
stars have steadily increasing periods. This is illustrated in Fig. 4 where we show a histogram of 
the (3 distribution for the stars that have a well defined O — C fit with monotonic changes exceeding 
the period noise. Of the 71 RRab stars listed in Table 4 and Table 5, 27 were excluded from the 
plot. Four (V9, 11, 56, 59) have ambiguous O — C curves (all of them exhibit Blazhko behavior) 
while for eight others (V91, 106, 132, 134,139, 144, 149, 150), there are not enough data to draw any 
firm conclusion about their long term period changes. In addition, those stars that did not show 
significantly monotonic changes in their periods according to the criterion of Eq. 5 are also omitted. 
V56 and V84 were also excluded because their radial velocities indicate that they are non-members 
(Liller and Tokarz 1981; van Leeuwen et al. 2000). 

Altogether, there remain 45 RRab stars for which monotonic period change rates seem to have 
been followed. Nine have decreasing periods; one of them (V104) has a steady period decrease at 
a very high rate (5.4 x 10~ 8 days/day). The period of this star, however, is unusually long for an 
ordinary RRab star (0.87 d). Thus its inclusion in the sample may be questioned. The other eight 
stars have periods decreasing at 10~ 10 and — 10~ 9 days/day magnitude rates. Period increases at 
similar rates have been detected in 30 of the 45 stars. For the other six, no change or very slight 
positive period changes were detected (0.0 < (3 < 10~ 10 days/day). The mean /3 of the 44 stars 
(omitting V104) is 4.1 x 10 -10 days/day (0.15 days/Myr), the median is 3.8 or 0.14 in the respective 
units. 

It was shown by Jurcsik (1998b), that in contrast to the large chemical inhomogeneity of 
the oo Centauri stars, the majority of the RRab stars have surprisingly homogeneous metallicity. 
Using the empirical formula that relates [Fe/H] to Fourier coefficients (Jurcsik and Kovacs 1996), 
metallicities of 48 RRab stars observed by Kaluzny et al. (1997) were derived. As no data were 
given in Jurcsik (1998b) we list the physical parameters of RRab stars calculated from the light- 
curves (Jurcsik 1998a; Kovacs and Jurcsik 1997, 1996; Jurcsik and Kovacs 1996) in Table 7 (see 
also Sect. 4). In this sample, 40 stars have similar metallicities; their [Fe/H]= —1.54 ±0.08 dex. Of 
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these 40, reliable O — C could be determined for 39 stars. Considering the period changes of only 
these variables, the tendency for the periods to increase is also clear. Positive (3 values have been 
detected in 24 stars, only 4 have negative ones, 1 does not show any period change and irregular 
period changes are observed in 10 cases. There was one star (V99) found to be more metal poor; 
its steady period increase rate is j3 = 46.0 x 10~ 10 days/day. Of the 7 relatively metal rich stars 
excluding the non-member V84 and the peculiar V104, 3 have shown steady period increases while 
only one is with decreasing period and the period change of one star is not monotonic. 

Regarding the Blazhko type behavior, we have found that 14 out of the 71 RRab stars clearly 
exhibit some type of modulation in their light-curves. Excluding the non-member V56, there is 
no sign of monotonic period change in five of the stars, while out of the eight variables showing 
monotonic period variation, a decreasing period has been observed for only one (V120). Although 
the periods of the w Centauri RRab stars span a very wide range (0.47 — 0.89 d), all the Blazhko 
variables have periods less than 0.7 d. 

In Fig. 5 the distribution of the normalized period changes (a) is plotted as a function of the 
period. Stars without a definite sign for their period changes are arbitrarily placed at a = —70. 
These stars are also found only among the shorter period variables, most of them have periods within 
a very narrow (0.59 — 0.67 d) period range. This unusual coupling of the periods of the Blazhko 
variables with the periods of stars showing random period variations might probably indicate some 
kind of physical relation between these phenomena. 

In conclusion, we have found that, although random variations dominate the period changes 
of 14 stars, period increases dominate in most of the remaining stars. Among the RRab variables 
exhibiting steady period changes, only 9 have decreasing periods while 34 have increased their 
periods during the past century. 

3.2. RRc stars 

We have constructed O — C plots for 48 RRc type variables, but out of which 26 must be taken 
with caution because of cycle number ambiguity or not enough time coverage. The RRc sample 
shown in Figs. 2a-h clearly indicates that the period changes of RRc stars are more complex than 
those of the fundamental mode (RRab) stars. Out of the RRc stars with well defined O — C curves, 
period increases occur in 10, period decreases in 4, random period changes in 7 cases, and 1 star 
does not show any period change. The equivalent numbers for the 26 stars with less certain O — C 
-s are 4, 3, 17, and 1, respectively, if V168 is not included because it is not a radial velocity member. 

A thorough examination of the light-curves of the RRc stars revealed the existence of some 
kind of modulation in at least 8 out of the 48 first overtone variables. This is a similar percentage 
(17%) to the occurrence of Blazhko behavior in the case of RRab stars (20%). 

We call attention to a special group of variables V47, V68, and V123 (showing also light- 
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curve variation) which we included in the RRc group because of their amplitudes and light-curves 
shapes, although it was shown in Clement and Rowe (2000) that their (f2i Fourier parameters are 
anomalously large. Geyer and Szeidl (1970) also noticed that some of the RRc stars define a bluer 
ridge in the B— Y.vs.logP plane, and among others V47 and V123 were found to belong to this 
group. 

The periods of these three stars are unusually long (0.47 — 0.53 d), which would rather indicate 
fundamental mode oscillation and their mean magnitudes are among the brightest of the RRc 
sample. V47 and V68 are 0.1 and 0.2 mag brighter than any of the others and the mean magnitude 
of V123 corresponds to that of the brightest among the other variables. Though the period changes 
of these stars are irregular according their uap /\(3\ parameter, in V68 and V123 significant period 

At 

increases (1.9 and 5.5 d/Myr, respectively) have occurred during the past century. Thus we may 
be seeing an indication of their more evolved status, although this is not as marked as in the case 
of the RR Lyrae-like stars discussed in the next section (3.3). These latter stars are typically more 
than 0.5 mag brighter than the RR Lyrae stars. There is, however a discernible similarity between 
the light-curves of these stars and that of the faintest evolved star (V92) even though it has a much 
longer (1.34 d) period. At present we can safely conclude that further studies, especially accurate 
color determinations are needed in order to clarify the true nature of these peculiar objects. 

Omitting these 3 variables and stars with too short O — C segments, a similar plot of the a 
distribution of the RRc variables as seen in Fig. 5 for the RRab stars is shown in Fig. 6. Again, 
stars with irregular period changes tend to be found among the longer period RRc stars where the 
modulated (Blazhko) variables are. 

3.3. RR Lyrae-like stars 

Wehlau and Bohlender (1982) made a thorough analysis of period changes of BL Her stars in 
globular clusters. They classified 5 of the u> Centauri variables as BL Her type and found period 
increase rates in the range of 10 -9 — 10 -8 days/day for all of them. 

The name BL Her is no longer used for these longer period RR Lyrae-like stars (short period 
population II Cepheids), but the classification scheme has not yet been settled (see e.g. Sandage 
et al. (1994); Nemec et al. (1994)). In this category, we consider the AHB1 stars (according to 
the nomenclature of Gautschy and Saio (1996)). AHB1 stars are 'above horizontal branch' stars 
which are either significantly brighter (0.5 — 1 mag) or have longer periods than the cluster RR 
Lyrae stars. Thus we also add V52 to this group because it is 0.5 mag brighter than other RRab 
stars with similar period (0.66 d). Based on its period we may also include V104 (P=0.87 d, 
(3 = —5.4 x 10~ 8 days/day) into this group, however its very large monotonic period decrease rate 
represents a disturbing extreme in any of the groups. 

We also included V43 and V48 in our study, although no new observations were available. 
The period change rates we have determined for these stars agree within the errors to those of 
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Wehlau and Bohlender (1982), showing the compatibility of our method with other period changes 
studies. The addition of new (Kaluzny et al. 1997) data to the earlier measurements strengthens 
the dominance of the monotonic period increases for V60,V61 and V92. Our solutions for V60 and 
V92 are, however, somewhat different from the earlier results (Wehlau and Bohlender 1982) mainly 
because we also take direct period determinations into account and this yields even larger period 
increase rates. 

Thus the period increases observed in these stars are in agreement with their evolved HB 
status. In the last phase of the HB, during evolution towards the asymptotic giant branch, large 
period increases are expected. 

4. Conclusion: Comparison with evolutionary model predictions 

We have shown in Sect. 3.1 that most of the RRab stars have increased their periods by some 
10~ n — 10~ 9 days/day rate during the past century. This is in agreement with evolutionary model 
predictions for the late, redward phase of the horizontal branch evolution. The period changes 
along the post zero age horizontal branch, as derived from the Dorman (1992) models and from 
the linear pulsation period formulae of Kovacs and Buchler (1994), range typically from — 1CT 10 
to +10 -9 days/day. Within the instability strip there is no metallicity and mass value for which 
the evolutionary period decrease rate is larger than some 10~ n days/day, but period increases of 
the order of 10 -9 days/day do occur in the redward phase of the blue loops in the HB evolution, 
especially for the lower mass stars. 

The physical parameters derived from the light-curves (Jurcsik 1998b, and Table 7) indicate 
that most of the RRab stars in to Centauri belong to a very homogeneous population, and are 
already in the evolved (redward) phase of their HB evolution. Thus, their comparison with model 
predictions of unique chemical composition models give adequate results, in spite of the fact that 
large chemical inhomogeneities have been detected in u) Centauri stars. 

In Fig. 7 horizontal branch evolutionary tracks are shown for oxygen enhanced models with 
[Fe/H]= -1.48 dex (Dorman 1992). Tracks for four different masses (0.66, 0.64, 0.62, O.6OM ) are 
shown. The position of the 40 RRab stars which were found to have [Fe/H]= —1.54 ± 0.08 dex 
(Jurcsik 1998b) are also given. The luminosities and temperatures of these variables (given in 
Table 7) were determined by using the empirical relations between the physical parameters and the 
Fourier coefficients of the light curves (see Jurcsik (1998a) and references therein). The empirical 
logL and logT scales are shifted by 0.10 and 0.016, respectively, in order to reach agreement with 
evolutionary results (Jurcsik and Kovacs 1999). This transformation yields a mean pulsational 
mass of 0.64M© with 0.02M© rms scatter for these 40 variables. 

Linear pulsational periods have been calculated for different locations along the tracks where 
the fundamental mode period falls within the 0.4 — 1.0 d range. Then, evolutionary period changes 
have been calculated using these periods and the evolutionary time scales of the Dorman (1992) 
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models. They are given in units of 1CT 10 days/day in Fig. 7. The distribution of the observed j3 
values shown in Fig. 4 for those RRab stars which have well defined monotonic period change rates 
is in very good agreement with the evolutionary (3 values calculated within the region occupied by 
these stars. 

Summarizing our results, we have found that, for those uj Centauri RRab stars that exhibit 
dominantly monotonic period changes, the period change rates can be brought into very good agree- 
ment with horizontal branch evolutionary model predictions. Thus the observed period changes of 
these stars may directly indicate evolutionary effects. However, irregular period changes and one 
very large period decrease rate (similar to those observed in other globular clusters) also occur in 
uj Centauri and no satisfactory explanation for these phenomena has yet been given. The much 
less regular period change behaviour of the first overtone RRc stars is also not yet understood. 

We are grateful to the referee, Emilia Belserene, for useful comments and suggestions. We 
also thank E. Belserene for sending some of her unpublished magnitudes. C. Clement would like 
to thank Peter Hanson and Tom Wells for making the iris photometer measurements. Fruitful 
discussions with B. Szeidl are highly appreciated. We are also grateful for financial support for this 
work which was partially provided by the Hungarian OTKA grants T-024022 and T-030954 and 
by the Natural Sciences and Engineering Research Council of Canada. 
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Fig. 1. — Results for RRab stars. O — C (left panels), directly measured period values (center) 
and folded light-curves of all the photometric data (right), after the times of the observations were 
transformed according to the polynomial fit of the O — C curve. 2 400 000 has been substructed 
from the JDs. Dashed lines in the O — C plots indicate an uncertainty in cycle number counts 
or in the shape of the O — C fit. The filled circles denote O — C points derived from the original 
photometric data and the open circles denote O — C values adopted, or estimated from folded light 
curves. The notation 'BP indicates evidence of light curve modulation. The period fits in the center 
panels correspond to the derivative of the O — C fit and not to actual fits to the periods. The filled 
circles represent direct period determinations. The magnitude scale in the right panel is CCD V if 
the star was measured by Kaluzny et al. (1997) and photographic m pg of Martin (1938) for all the 
other stars. 

Fig. 2. — Results for RRc stars. See details in caption Fig. la. 

Fig. 3. — Results for RR Lyrae like stars. See details in caption Fig. la. 

Fig. 4. — Histogram of period change rates (/?) of the 44 RRab stars that have unambiguous O — C 
curves and monotonic period changes dominating over irregular period noise. The bulk of these 
RRab stars have continuously increased their periods by some 10~ 10 days/day rate during the 
century long time-base of the observations. 

Fig. 5. — Normalized period change rates (a) vs. periods of RRab stars. Open circles denote 
Blazhko variables. V104 (a = -625 x 10" 10 /day), the non- members (V56, V84), and stars that do 
not have long enough continuous coverage on their O — C curves are not shown. Those variables 
which have some ambiguity in their O — C fits are included because their dominant (random or 
monotonic) changes can be safely determined. Alpha values of those variables which exhibit random 
period changes without significant monotonic contribution are arbitrarily set to —70 just to show 
their distribution with period. All of these stars and the ones with decreasing periods belong to 
the shorter period group of the RRab sample, similar to the Blazhko type variables. 

Fig. 6. — The same plot as Fig. 5 for the RRc stars. The anomalously long period stars (V47,68,123) 
and V151,V163 and V168 are not shown. Stars with random period changes are arbitrarily set to 
a = -150. 

Fig. 7. — Horizontal branch evolutionary tracks of the oxygen enhanced, [Fe/H]= —1.48, M = 
0.60, 0.62, 0.64, O.66M models within the region of the instability strip (Dorman 1992). Evolu- 
tionary period changes along the tracks ({3 in 10~ 10 days/day unit) are given. Filled circles denote 
the 40 variables with [Fe/H]= —1.54 ± 0.08 dex and M = 0.64 ± O.O2M pulsational mass values. 
Long term period changes for 39 of these stars have been determined: period increase with (3 = 6.6 
mean rate for 24 stars, decreasing periods for 4 stars, and random period changes for 10 stars. One 
star does not show any period change. The largest period increase rates (j3 > 20.0) are observed in 
the two most luminous stars. The mean period increase rates for the stars with positive (3 values 
are in very good agreement with the evolutionary results. 



Table 1. Yale-Columbia observations of lo Centauri variables. The complete version of this table 
is in the electronic edition of the Journal. The printed edition contains only a sample. 



J.D.-2 400 000 


V3 


V7 


V8 


V10 


V12 


V13 


V14 


V15 


V16 


V18 


V19 


32234.585 


14.32 


14.42 


14.00 


14.86 


14.70 


14.53 


14.50 


14.15 


15.05 


15.36 


15.10 


32266.500 


14.18 


14.52 


14.68 


14.98 


14.58 


14.56 


14.99 


14.72 


14.64 


14.48 


14.72 


32278.462 


14.57 


15.03 




14.92 


14.49 


15.22 


14.59 


14.36 


14.98 


15.02 


14.77 


32295.414 


14.88 




15.41 


14.93 


14.48 


14.38 


14.48 


14.22 


14.54 


15.29 


15.24 


32296.410 


15.01 


15.10 


15.42 


14.82 


14.96 


15.09 




14.67 




14.83 


15.06 



Table 2. La Campanas observations of to Centauri variables. The complete version of this table 
is in the electronic edition of the Journal. The printed edition contains only a sample. 



J.D.-2 400 000 


V3 


V4 


V5 


V8 


V9 


V10 


V13 


V16 


V18 


V19 


V20 


V22 


V23 


V24 


41446.739 


15.08 


15.29 


15.32 


15.32 


14.35 


14.96 


14.14 


14.99 


14.79 


15.32 


14.85 


15.01 


15.32 


14.99 


41447.555 


15.26 


14.34 


14.54 


14.16 


15.34 


14.97 


14.80 


14.41 


15.29 


14.76 


15.17 


15.14 


14.30 


14.66 


41447.655 


14.38 


14.92 


15.13 


15.00 


15.38 


14.51 


15.07 


14.61 


15.25 


15.28 


15.28 


14.69 


15.02 


15.10 


41448.507 


14.26 


15.22 


14.26 


15.33 


15.29 


14.46 


15.11 


14.60 


14.30 


15.15 


14.34 


14.41 


15.04 


14.67 


41449.475 


14.25 


14.48 


15.41 


15.44 


15.10 


14.80 


14.39 


14.96 


15.26 


15.41 


15.20 


15.05 


15.35 


14.92 



Table 3. ADH Schmidt observations of uj Centauri variables. The complete version of this table 
is in the electronic edition of the Journal. The printed edition contains only a sample. 



J.D.-2 400 000 


V4 


V10 


V14 


V18 


V19 


V22 


V24 


V36 


V56 


37792.289 


14.85 


14.60 


14.81 


15.26 


14.74 






14.99 


15.27 


37792.315 


14.75 


14.35 


14.74 


15.05 


14.84 


14.57 


14.25 


14.89 


15.35 


37792.358 


14.82 


14.25 


14.96 




14.45 


14.31 


14.32 


14.99 


15.01 


37793.341 


13.99 


14.74 


14.50 


15.11 


14.73 


14.93 


14.57 


14.59 


15.35 


37793.368 


14.17 


14.45 


14.52 


15.16 


15.00 


14.88 


14.77 


14.62 


15.49 



Table 4. Photometric data of RR Lyrae and evolved 
variables in to Centauri. 



Var. a References' 1 



SH 


K 


1 


2 


3 


4 


5 


6 


7 


8 


9 


type 


3 


184 


X 


X 


X 





xo 


— 


X 


X 


X 


ab 


4 


99 


X 


X 


X 


X 


— 


— 


X 


X 


X 


ab 


5 


101 


X 


X 


X 





— 


— 


X 


X 


— 


ab 


7 


87 


X 


X 


X 





xo 


X 


— 


X 


— 


ab 


8 


199 


X 


X 


X 





xo 


— 


X 


X 


— 


ab 


9 


183 


X 


X 


X 


— 


— 


— 


X 


X 


X 


ab 


10 


98 


X 


X 


X 


X 


X 


— 


X 


X 


X 


c 


11 


126 


X 


X 


— 


c 


— 


— 


— 


X 


X 


ab 


12 


125 


X 


X 


— 


— 


X 


— 


— 


X 


X 


c 


13 


188 


X 


X 


X 





X 


— 


X 


X 


X 


ab 


14 




X 


X 


X 


X 


X 


— 


— 


— 


— 


c 


15 


124 


X 


X 


— 





xo 


— 


— 


X 


X 


ab 


16 


205 


X 


X 


X 





X 


— 


X 


X 


— 


c 


18 


200 


X 


X 


X 


X 


xo 


— 


X 


X 


— 


ab 


19 


192 


X 


X 


X 


X 


X 


— 


X 


X 


— 


c 


20 


163 


X 


X 


X 


— 


X 


— 


X 


X 


X 


ab 


21 


97 


X 


X 


— 


— 


X 


— 


— 


X 


X 


c 


22 


190 


X 


X 


X 


X 


— 


— 


X 


X 


X 


c 


23 


146 


X 


X 


— 


— 


xo 


— 


X 


X 


X 


ab 


24 


189 


X 


X 


X 


X 


X 


— 


X 


X 


— 


c 


25 


113 


X 


X 


— 


— 


xo 


— 


X 


X 


— 


ab 


26 


120 


X 


X 




— 


X 


— 


X 


X 


— 


ab 


27 


114 


X 


X 


— 


— 


X 


— 


X 


X 


X 


ab 


30 


92 


X 


X 












X 


X 


c 


32 


82 


X 


X 


— 


— 


xo 


— 


X 


X 


— 


ab 


33 


182 


X 


X 


X 





xo 


— 


— 


X 


— 


ab 


34 


90 


X 


X 


X 





X 


— 


X 


X 


X 


ab 


35 


78 


X 


X 


X 


— 


X 


— 


X 


X 


X 


c 


36 


86 


X 


X 


X 


X 


— 


— 


— 


X 


— 


c 


38 


170 


X 


X 


X 





xo 






X 




ab 


39 


197 


X 


X 


X 





X 


— 


X 


X 


X 


c 


40 


107 


X 


X 




— 


xo 


— 


X 


X 


X 


ab 


41 


129 


X 


X 




— 


X 


— 


X 


X 


X 


ab 


43 




X 


X 


X 








X 






evolved 


44 


88 


X 


X 


X 





xo 




X 


X 


X 


ab 


45 


179 


X 


X 


X 







X 


X 


X 




ab 


46 


180 


X 


X 


X 





X 


X 


X 


X 




ab 


47 


185 


X 


X 




c 








X 


X 


c 


48 




X 


X 


X 








X 






evolved 


49 


165 


X 


X 


X 




X 




X 


X 




ab 


50 


187 


X 


X 


X 





X 




X 


X 




c 


51 


147 


X 


X 






X 




X 


X 


X 


ab 


52 


116 


X 


X 












X 




evolved 


54 


74 


X 


X 







xo 




X 


X 




ab 


55 




X 


X 







xo 


X 








ab 


56 




X 


X 


X 


X 


X 




X 






ab 


57 


206 


X 


X 


X 





xo 




X 


X 




ab 


58 


100 


X 


X 




X 


X 






X 


X 


c 


59 


93 


X 


X 












X 


X 


ab 


60 


103 


X 


X 


X 








X 


X 




evolved 


61 


164 


X 


X 


X 








X 


X 




evolved 



Table 4 — Continued 



SHK 123456789 type 

62 117 X X - - X - - X X ab 

63 X X X XO Y - - - ab 

64 XXXOX — X — — c 

66 76 XX-----XX c 

67 77 XX — X — — XX— ab 

68 73 XX — XX — — X— c 

69 XX-O-X--- ab 

70 151 XXXXX-XXX c 

71 157 XX-----XX c 

72 XXXOX---- c 

73 XX-O-X--- ab 

74 85 X X X XO - X X - ab 

75 84 XX-OX--X- c 

76 79 XXXOX-XX- c 

77 80 XX — C — — — XX c 
79 204 X X X X - - X - ab 

81 195 XXXOX-XX- c 

82 193 XXX-X-XX- c 

83 83 XXXXX-XX- c 

84 174 X X - X X - X - ab 

85 176 X X - X Y - X - ab 

86 158 X X - - X - X X X ab 

87 144 XX-----XX c 

88 210 X X - - XO - - X - ab 

89 143 XX-----XX c 

90 141 X X ----- X - ab 

91 XX------- ab 

92 72 XXXO-XXX - evolved 

94 XX-OX---- c 

95 178 XXXOXXXX- c 

96 XX----X-- ab 

97 159 X X - - X - - X X ab 

98 160 XX-----X- c 

99 137 X X - - X - - X X ab 

100 X X - - X - - - - ab 

101 194 X X X - X - - X - c 

102 153 X X - - X - - X X ab 

103 XX--X---- c 

104 201 X X X XO - X X X ab 

105 198 X X X X X - X X - c 

106 X X ------ - ab 

107 154 X X - - X - - X X ab 

108 211 X X - - XO - - X - ab 

109 133 X X - - X - X X - ab 

111 136 X X ----- X X ab 

112 150 X X ----- X - ab 

113 128 X X - - X - X X - ab 

115 89 X X X - - - X X X ab 

116 127 X X ----- X - ab 

117 94 XX — — X — — XX c 

118 109 X X ----- X - ab 



Table 4 — Continued 



Var. a References 1 " 



SH 


K 


1 


2 


3 


4 


5 


6 


7 


8 


9 


type 


119 


106 


X 


X 












X 


X 


c 




1 CiA 


V 


y 
. \ 
















dD 


121 


105 


X 


X 






X 






X 




C 


122 


102 


X 


X 






X 




X 


X 




ab 


123 


169 


X 


X 







X 






X 


X 


c 


124 


168 


X 


X 




X 


X 






X 




c 


125 


166 


X 


X 


X 





X 






X 




ab 


126 


207 


X 


X 




X 


X 






X 




c 


127 


175 


X 


X 


X 





X 


X 


X 


X 




c 


128 


91 


X 


X 






xo 






X 




ab 


130 


75 


X 


X 




c 








X 




ab 


131 


110 


X 


X 












X 


X 


c 


132 




X 


X 
















ab 


134 






X 







X 










ab 


139 


118 




X 












X 




ab 


144 


112 




X 












X 




ab 


149 






X 







X 


X 








ab 


150 






X 






X 










ab 


151 






X 






X 










c 


155 


145 




X 






X 






X 


X 


c 


160 






X 







X 


X 








c 


163 






X 







X 










c 


168 


181 










X 






X 




c 



a SH refers to Sawyer Hogg (1973) and Clement (2000) catalogue numbers, 
K to the Kaluzny et al. (1997) IDs. 

b X: photometry; Y: photometry too few to define O — C value; 0: folded 
light curve or O — C data available; C: only folded light curve available without 
parameters allowing to determine O — C . 

References. — (1) Bailey (1902); (2) Martin (1938); (3) Dickens and Saun- 
ders (1965) or Dickens and Carey (1967) for evolved stars; (4) Geyer and Szcidl 
(1970) data published in this paper and read from their folded light-curves; (5) 
Belserene (1964) data published in this paper and Belserene (1964) O — C 
values for JD 2423530, 2437454 and 2437727; (6) Sturch (1976, 1978); Eggen 
(1961) for V92; (7) C. Clement observations published in this paper or Wehlau 
and Bohlcnder (1982) for evolved stars; (8)-(9) observations of different OGLE 
fields Kaluzny et al. (1997). 



Table 5. O — C and period data of uj Centauri variables. The complete version of this table is in 
the electronic edition of the Journal. The printed edition contains only a sample. 



Var. P a (d) a order of the polynomial fit to the O — C 

mean JD b O - C (d) first JD b < c last JD b - c N c - d Am c - C Rcf. f 

mean JD b measured period period range 



V003 0.84f22900 6 



f3f53.06 


-0.036 


12234.7 


13411.6 


43 


0.88 


1 


f4f74.72 


-0.064 


13652.6 


14518.5 


61 


0.88 


1 


23530.00 


-0.383 


0.0 


0.0 





0.00 


5 


iU'ifJ l.O I 


-0.449 






263 


-0.24 


2 


27907.03 


-0.453 


27891.4 


27988.0 


135 


-0.24 


2 


32463.33 


-0.486 


32234.6 


32736.2 


24 


-0.28 


5 


35670.10 


-0.503 


35607.0 


35682.0 


27 


-0.28 


5 


36269.99 


-0.539 


36259.2 


36292.2 


13 


-0.28 


5 


37298.90 


-0.498 


36754.2 


37491.2 


36 


-0.36 


3 


37454.00 


-0.481 


0.0 


0.0 





0.00 


5 


37811.00 


-0.485 


0.0 


0.0 





0.00 


4 


38084.97 


-0.457 


38062.1 


38118.1 


12 


-0.28 


5 


38885.84 


-0.466 


38876.1 


38909.0 


11 


-0.28 


5 


39629.88 


-0.453 


39563.1 


39684.9 


28 


-0.28 


5 


41669.38 


-0.385 


41446.7 


41874.5 


43 


-0.34 


7 


42359.68 


-0.355 


42074.8 


42594.5 


53 


-0.34 


7 


43520.57 


-0.303 


43273.5 


43699.7 


73 


-0.34 


7 


44194.27 


-0.272 


44016.6 


44378.8 


53 


-0.34 


7 


44993.61 


-0.237 


44750.6 


45470.5 


79 


-0.34 


7 


46244.63 


-0.175 


45786.8 


46562.7 


54 


-0.34 


7 


49524.49 


-0.016 


49515.5 


49525.7 


28 


-0.02 


8 


49831.88 


0.000 


49818.5 


49922.5 


201 


0.00 


9 


13752 


0.84120425 


0.84118370 


0.84122730 








26945 


0.84122507 


0.84121650 


0.84123400 








36544 


0.84122675 


0.84121620 


0.84124020 








39105 


0.84125025 


0.84123810 


0.84126290 








42685 


0.84126670 


0.84125840 


0.84127480 








45129 


0.84126682 


0.84125890 


0.84127450 








49794 


0.84127000 


0.84125740 


0.84128260 









a The O — C values are calculated according to this period. 
b JD-2 400 000 

c 0.0 values are for references when the original photometric data are not available. 
d Number of data points used. 

e Magnitude-shift applied to vertically match data to the normal curve. 
f Reference numbers are the same as in Table 4. 



Table 6. Period changes parameters of u Centauri variables. 



Var. a P a (d) Bl b 0-C c [Fc/H] d f3 =< ^ > a = /3/P a ^ v±r/\/3\ 

mean period d/10 10 d d/Myr l/10 10 d 



RRab stars 



3 


0.8412290 


1 


1 


1 


21.750 


0.794 


25.855 


1.975 


0.09 


4 


0.6273165 


1 


1 


1 


3.733 


0.136 


5.951 


0.649 


0.17 


5 


0.5152828 





1 




1.617 


0.059 


3.137 


0.832 


0.51 


7 


0.7130227 


1 


1 


1 


6.833 


0.249 


9.584 


0.901 


0.13 


8 


0.5212930 


1 


1 


1 


16.950 


0.619 


32.515 


0.696 


0.04 


9 


0.5233400 










32.314 


1.179 


61.745 


20.363 


0.63 


11 


0.5648060 










5.042 


0.184 


8.926 


0.008 


0.00 


13 


0.6690535 


1 


1 


1 


2.528 


0.092 


3.778 


2.725 


1.08 


15 


0.8106220 


1 


1 


1 


20.406 


0.745 


25.173 


0.004 


0.00 


18 


0.6216686 


1 


1 


1 


1.217 


0.044 


1.957 


1.585 


1.30 


20 


0.6155540 


1 
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0.15 
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-625.112 


1.460 
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a *: not radial velocity members (Liller and Tokarz 1981; van Leeuwen et al. 2000). 

b 0: light-curve showing any type of modulation. 

c 0: uncertain O — C fit; not enough data to draw firm conclusion about long term period changes. 

d 1: [Fc/H]= —1.54 ± 0.08 dex; 0: less metal poor; 2: more metal poor stars (see Table 7). 



Table 7. Physical parameters of RRab stars determined from their light-curves. 
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Note. — The physical parameters are calculated using the formulae 
given in Jurcsik (1998a), the LogL and LogT scales are shifted by 0.1 
and 0.016, respectively, in order to reach agreement with the Dorman 
(1992) evolutionary models (Jurcsik and Kovacs 1999). 
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